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The Cucchieri-Mendes bounds for the gluon propagator are discussed for the four dimensional 
pure gauge SU(3) theory. Assuming a pure power law dependence on the inverse of the lattice 
volume, the lattice data gives a vanishing zero momentum gluon propagator in the infinite volume 
limit. Our investigation shows that this result is robust against variations of the lattice volumes and 
corrections to the power law. Moreover, it is consistent with the Gribov-Zwanziger horizon condition 
and contradicts the SU(2) analysis, which assumed the same type of dependence on the inverse of 
the volume. Our analysis considers also more general ansatze that, although not conclusive, open 
the possibility of having -D(O) 7^ in the infinite volume limit. A solution to this puzzle requires 
further investigations. 



I. INTRODUCTION AND MOTIVATION 

Gluon confinement criteria [l[ make precise state- 
ments about the behavior of the zero momentum Landau 
gauge gluon propagator I?(0). In particular, the Gribov- 
Zwanziger horizon condition predicts a van- 
ishing propagator at zero momentum. Moreover, in it 
was shown that, in the minimal Landau gauge, the lattice 
zero momentum gluon propagator vanishes in the infinite 
volume limit and that same prediction holds for all zero 
momentum gluon correlation functions. For the gluon 
propagator, these predictions are not in line with early 
calculations using Dyson-Schwinger equations which re- 
quired an infinite zero momentum gluon propagator to 
explain quark confinement 0, i.e. a linearly rising 
interquark potential. 

First principles calculations of the gluon propagator 
rely either on Schwingcr-Dyson equations or on lattice 
QCD simulations. 

In what concerns the Schwinger-Dyson equations, re- 
cent solutions follow in two distinct classes: i) solutions 
which have a vanishing propagator at q = (for a recent 
review see [loj): ii) solutions with a finite but non van- 
ishing value for the propagator at zero momentum 
The reader should be aware that the Schwinger-Dyson 
equations are an infinite tower of equations relating the 
Green's functions of QCD. Finding a solution of the in- 
finite set of equations requires truncating the tower of 
equations and parametrizing some of the vertices. 

On the other hand, four dimensional lattice simula- 
tions for the gluon propagator have been done with SU(2) 
and SU(3) gauge theories (see, for example, El, 
EHIlEHliS,!!!,!!!,!!!] and references therein). All 
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lattice simulations show a finite -D(O) 7^ 0. Furthermore, 
recent simulations using huge lattice volumes (27 fm)'^ for 
SU(2) [13 and (13.2fm)4 for SU(3) [13 give a propaga- 
tor which goes to a constant in the infrared limit. More, 
the raw lattice data shows no suppression of D[q^) in 
the infrared limit. However, the lattice data shows a de- 
creasing D{Q) as the lattice volume increases and naive 
extrapolations of I?(0) to the infinite volume, which typ- 
ically assume a linear behaviour, always give a finite non 
vanishing value around 6—10 GeV~^. These results are 
in contradiction with the prediction of a vanishing D{Q) 
in the infinite volume limit 

This contradiction can be viewed in two ways, either 
the simulated volumes are still too small or the finite vol- 
ume effects provide sizable corrections to the gluon prop- 
agator in the infrared limit which are not under control. 
Giving the typical hadronic scales ~ 1 fm, the scale for 
refiection positivity violation ~ 1 fm, and the volumes 
of the recent lattice simulations ^ 13 fm or larger, it is 
harder to believe that such volumes are not sufficiently 
large. 

To overcome the problem of the finite size effects, in 
[ist a method was suggested which suppresses the finite 
volume effects. Indeed, the simulations reported in the 
article use a set of different large asymmetric lattices 
and for all the lattices investigated in [ll|, D{0) is fi- 
nite and non vanishing. The reader should be aware that 
despite the observed differences on the gluon propagator, 
clearly due to finite volume and asymmetry effects p^ . 
the method discussed in [ll| to access the infrared prop- 
erties of the gluon propagator give compatible results for 
all lattices investigated in the article. Furthermore, in 
[ist no a priori assumptions about the value of -D(O) in 
the infinite volume limit were made. The new method 
discussed in that work provides consistent results for all 
lattices and points always to a vanishing D(Q) in the in- 
finite volume limit. 
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In a previous article 



the same authors arrived 
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at a similar conclusion based on extrapolations towards 
the infinite volume of the fitted propagators. Presently, 
arc the only lattice simulations which suggest a 
result in agreement with the Gribov-Zwanzigcr horizon 
condition and with the results of [Q] for the giuon propa- 
gator but they required a closer look at the finite volume 
effects. 

In [2^, the authors derived rigorous upper and lower 
bounds for the zero momentum gluon propagator of lat- 
tice Yang-Mills theories in terms of the average value 
of the gluon field. The bounds follow directly from the 
Monte Carlo approach to lattice simulations. More in- 
teresting is the scaling analysis which the authors have 
performed for the SU(2) Yang-Mills theory. The scaling 
analysis allows to conclude in favor or against a 13(0) = 
in the finite volume. For SU(2) the results show a D{0) 
going to zero for the two dimensional theory, in agree- 
ment with a previous investigation p^l, but a non van- 
ishing D{0) for three and four dimensional formulations. 

The aim of this article is to perform a similar investi- 
gation for the four dimensional lattice SU(3) Yang-Mills 
theory and, hopefully, to check the infinite volume value 
for D{0). To achieve our goal we use different sets of 
lattices. The scaling analysis performed to each set of 
lattices concludes always in favor of a D{Q) when 
V — > +00, in agreement with Gribov-Zwanziger scenario 
and [i,[l3,[23. Note that, although in [25i,l^ the SU(2) 
and SU(3) lattice gluon propagators seem to be equal up 
to momenta below 1 GeV, the results of scaling analy- 
sis for the two theories (see [2^) suggests that the vol- 
ume dependence is not the same for both theories. Our 
investigation docs not provide a final answer on these 
questions and further investigations are required to un- 
derstand if there is a difference between the two theories 
and the reason for such a difference. 

The paper is organized as follows. In section |TT] wc 
briefly review the Cucchieri-Mendes bounds and discuss 
the different sets of lattices used here. Finally, in section 
mil we show the results of the bounds and discuss the 
scaling analysis. 

II. DEFINITIONS AND LATTICE SETUP 

The Cucchieri-Mcndes bounds relate the gluon propa- 
gator at zero momentum -D(O) with 

where d is the number of space-time dimensions, Nc the 
number of colors, 

^m(O) = ^EA"(^) (2) 

X 

and A'^{x) is the a color component of the gluon field 
in the real space. According to [2^, the D{0) is related 



TABLE L Symmetric lattices (left). The first column shows 
the lattice side in fm, the second column reports the num- 
ber of lattice points, the third the number of configurations 
generated and the last column the first non- vanishing momen- 
tum. For the conversion to physical units we used = 1.94 
GeV [2^. Asymmetric lattices (right). The table shows the 
lattices and the number of configurations for each lattice. 



La 


L* 


# 


Qmin 


xT 


# 


(fm) 




Conf. 


(MeV) 




Conf. 


1.63 


16'* 


52 


757 


8^ X 256 


80 


2.04 


20* 


72 


607 


10^ X 256 


80 


2.45 


24* 


60 


506 


12^ X 256 


80 


2.86 


28* 


56 


434 


14^ X 256 


128 


3.26 


32* 


126 


380 


16^ X 256 


155 


4.90 


48* 


104 


254 


18''* X 256 


150 


6.53 


64* 


120 


190 







with M (0) by 

(A/(0))2 < ^ <d{N^-l){M{Of). (3) 

In the last equation ( ) means Monte Carlo average over 
gauge configurations. The bounds in equation Q are a 
direct result of the Monte Carlo approach. The interest 
on these bounds comes from allowing a scaling analy- 
sis which can help understanding the finite volume be- 
haviour of D{Q): assuming that each of the terms in in- 
equality ^ scales with the volume according to A/V"', 
the simplest possibility and the one considered in [23| , an 
a > 1 for (M(0)^) clearly indicates that D{0) ^ as the 
infinite volume is approached. 

In our lattice investigation we used the SU(3) pure 
gauge Wilson action at /? = 6.0 for a number of different 
sets of lattices - see table HI The gauge configurations 
were generated with version 6 of the MILC code [l^l . 

For the 48'^ and 64"' lattices, the configurations were 
gauge fixed using an overrelaxation algorithm. For the re- 
maining lattices, the gauge fixing procedure was a Fourier 
accelerated steepest descent. See [l^] for more details on 
the definitions of the gluon field and propagator. 



III. RESULTS 

In figures [T] and [2] wc show {M{0)f, D{0)/V and 
d(iVp — l)(A/(0)^) for the symmetric and asymmetric sets 
of lattices, respectively, together with the fits to A/V". 
As for the SU(2) case [2^, the bounds are verified with 
£1(0)/^ being closer to d{N^ - l)(A/(0))2. 

The infinite volume limit of D{0) can, in principle, be 
resolved studying the behaviour of the upper bound in 
([5]) with the volume. As stated previously, this was done 
assuming that each of the terms in ^ follows a A/V" 
law. For the symmetric lattices in table I, the fits give 
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Symmetric Lattices 

L"*. L=| 16,20,24,28,32,48,641 

- I I I I I - 




FIG. 1: Cucchieri-Mendes bounds for the symmetric set of 
lattices reported in table IH 
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FIG. 2: The same as figure [T] for the asymmetric lattices 
reported in table |T] 




a xl 

(A/(0)) 0.5289(29) 2.3 

D{0)/V 1.0587(53) 2.4 

(M(0)2) 1.0555(57) 2.4 

where xt ~ X^/d-o.f. Despite the slightly high value for 
the X /d-o.f., the results for D{0)/V and (M(0)2) show 
that, in the infinite volume limit, D{0) = 0. Moreover, 
they show that the approach to the infinite volume limit 
is rather slow Dv{0) - y-0 059(5)^ 

A careful looking at figure [1] shows that the -0(0) data 
from the 24* lattice is slightly above the fitting curve 
and is responsible for such high x^ /d-o.f. values. Indeed, 
removing the 24* data and repeating the analysis one 



gets 

a xl 

{M{0)) 0.5267(29) 1.0 

D{0)/V 1.0558(55) 0.8 

(M(0)2) 1.0526(59) 1.0 

i.e. the a's from the two fits arc compatible within one 
standard deviation, confirming, in this way, the results of 
the first fittings. The reader should note that in all cases, 
a from the fits to D{0)/V only becomes compatible with 
1, which would allow a finite D(0) in the limit of infinite 
volume, within 10 standard deviations. 

For the asymmetric lattices, the fits reproduce the 
same results 

a xl 
(M(0)) 0.5320(68) 0.8 
D{0)/V 1.104(13) 0.5 
(Af(0)2) 1.062(14) 0.7 

but with a x^ /d.o.f. ~ 1. Again, D{Q) approaches its 
infinite volume rather slowly Dyifi) ^ y-o io(i). Note 
that the a's for the asymmetric set are not compatible, 
within one standard deviation, with the values computed 
using the symmetric set. However, to the best of our 
knownledge there is no reason why the scaling laws for 
symmetric and asymmetric lattices should be exactly the 
same. Nevertheless, they should provide the same infi- 
nite volume limit and this is indeed the case: both sets 
provide Doo(O) = 0. 

The reader should note that the a from the fits to 
D[Q)/V and (M(0)^) are above 1 by more than several 
standard deviations both for the symmetric and asym- 
metric sets of lattices. In this sense, these results show 
that D{0) = in the infinite volume limit [ssj . 

The observed slow approach of -D(O) to its infinite vol- 
ume limit could explain why the previous naive extrapo- 
lations to y — > +00 which assumed a linear dependence 
OTil/V work so well. 

The volumes considered in this work range from 
(1.63 fm)* up to (6.53 fm)^ are smaller than the vol- 
umes of the SU(2) analysis [2^, where the largest volume 
is ~ (27 fm)*. On the other hand, the SU(2) simulations 
reported have coarser lattices. For our SU(3) simulation 
a = 0.102 fm, to be compared with a ~ 0.21 fm for the 
SU(2) case; note that both studies use the Wilson action. 
How the finite volume and finite lattice spacing effects 
change the conclusions reported in the two simulations is 
an open question. 

Nevertheless, in order to try to have an estimate of 
the influence of the lattice volume on our results, in ta- 
ble |TT] we have performed the same fit to subsets of our 
symmetric lattices. The a values are, within one stan- 
dard deviation, stable and compatible with the previous 
results. In particular, one can observe that restricting 
oneself to volumes up to 28*, with a physical volume less 
than (3/m)*, one gets a values well above 0.5 or 1 - see 
also figure O 
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{Aim 


D{Q)/V 


d{N^, -l){M{Qf) 


Lattices 


a xl 


a 


xl 


a. 


xl 


16,20,28,32,48,64 


0.5267(29) 0.98 


1.0558(55) 


0.80 


1.0526(59) 


1.01 


20,28,32,48,64 


0.5261(38) 1.29 


1.0542(72) 


1.02 


1.0520(77) 


1.35 


28,32,48,64 


0.5216(52) 1.09 


1.049(10) 


1.27 


1.042(10) 


1.11 


32,48,64 


0.5178(60) 0.58 


1.043(12) 


1.41 


1.035(12) 


0.56 


16,20,28 


0.524(10) 0.14 


1.052(20) 


0.0055 


1.045(21) 


0.22 


16,20,28,32 


0.5330(62) 0.59 


1.063(12) 


0.24 


1.064(12) 


0.74 


16,20,28,32,48 


0.5312(41) 0.44 


1.0641(75) 


0.17 


1.0615(81) 


0.53 


16,20,28,32,48,64 


0.5267(29) 0.98 


1.0558(55) 


0.80 


1.0526(59) 


1.01 



TABLE IL Fitting subsets of the symmetric lattices to A/V"' , not taking into account the 24*. 



1.08 
1.07 
1.06 
B 1.05 
1.04 
1.03 
1.02 
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1 



<M(0)>" 

D(0)/V 

<M(Or> 
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1/V 



FIG. 3; a from fitting the lattice data after removing the 
largest lattices. On the x-axis, 1/V refers to largest lattice 
used in the fit, i.e. smallest 1/V means including largest 
lattices. For (A'/(0)), the figure shows 2a. 



One could also try to evaluate the contributions from 
the corrections to the leading behaviour. Having this in 
mind, the data was fitted to 



The fits to this functional form in (j4]) gives 



(4) 



(M(0)) 
D{Q)/V 

(Mm 



0.5232(53) 

1.050(10) 

1.047(11) 



1530 ± 1954 
2518 ±3808 
2689 ±4012 



1.1 
0.9 
1.2 



Note that these new set of a's are, within one standard 
deviation, compatible with the values computed assum- 
ing a A/V" behaviour, i.e. the analysis of the corrections 
to the leading scaling behaviour do not change the above 
results. Furthermore, w ~ within the statistical accu- 
racy of the simulation. The fits to the asymmetric lattices 
data, not reported here, lead to the same conclusions. 

The above analysis clearly supports a vanishing gluon 
propagator in the infinite volume, i.e. DadO) — 0. How- 
ever, as already discussed in [3 



the value of Doo(O) 



depends on the way we model the approach to the infi- 
nite volume. The reader should be aware that there is no 
theoretical guidance about the behaviour of these quan- 
tities with the volume. Up to this point, as in [23[, we 
have assumed the simplest behaviour, a power law in the 
volume. In order to test for the possibility of having a fi- 
nite and non- vanishing _Doc (0), one has to consider other 
fitting functions. In particular we tested the following 
one. 



C 
V 



A 

V"' 



(5) 



A finite and non- vanishing Doo(O) requires C ^ 0. For 
the full symmetric set of lattices considered in this work, 
we get fits with y^/d.o.f. > 2.5, i.e. the quality of the 
fits arc poor. However, if wc do not take into account the 
24^ data as done previously, the fits give 



A 

(M(0))2 208 ±228 
D{0)/V 263 ±247 
(il/(0)2) 6127 ±6645 



1.21(12) 
1.18(12) 
1.20(13) 



C 

33(7) 

48(14) 

1064(236) 



0.8 
0.7 
0.9 



(note that in the first line we are reporting (Af (0))^ and 
not (Af(0))). 

For the asymmetric set, we get 



(M(0))2 273 ±6020 
£1(0)/^ 539 ±6849 
(A/(0)2) 3957 ± 15273 



a C 

1.01(31) -190 ±6102 

1.02(35) -332 ±7151 

1.14(65) 788 ± 2891 



1.06 
0.67 
0.94 



So, while for the asymmetric lattices one gets C = 0, 
the analysis for the symmetric lattices allows for a finite 
non-zero Z?oo(0). In this sense, the results of ([5]) are not 
conclusive. 

As before, one could try to check for the stability of the 
fits considering subsets of symmetric lattices - see table 
mil The results shows that we only get C when our 
largest lattice 64^ is included in the analysis. All the 
results of fitting equation ^ suggest that either it does 
not make sense to use ([5]) or larger volumes are required. 
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We would like to recall that a similar analysis assuming 
A/V^ behaviour does not suffer from the same problem, 
i.e. the lattice data is well described by A/V" and the 
exponent is independent of the subset of lattices used in 
the analysis. 

For completeness, we have also considered a general- 
ization of ^ and fitted the lattice data to 



(6) 



(note that equation ([5]) corresponds to the particular 
case a = 1). A full analysis of this ansatz is not 
easy, given the non-linear nature of the fit. However, 
the task can be made easier fixing a to the values 
obtained from the fits to the pure power law and 
then evaluated /3. For the symmetric lattices, dis- 
carding the data from 24^, the results are as follows 



(M(0))2 93.6 ±1.4 
D{0)/V 152.2 ±2.1 
(A/(0)2) 83 ±46 



0.87 ± 1.90 
0.97 ±2.95 
0.79 ±1.89 



W2 xl 

240 ± 5061 1.24 

574 ±19718 1.01 

83 ± 1740 1.30 



Note that in all cases and within the present 
statistical accuracy, lu2 ~ 0. However, if 
we fix a = 0.9, which would imply an infi- 
nite -D(O), the x^/d.o.f. are again around one. 



(A//(0))2 3.11 ±1.6 
£'(0)/^" 4.2 ±2.6 
(M(0)2} 98 ±54 



/3 

0.216(46) 
0.207(43) 
0.212(46) 



^2 XI 

50.8 ±8.1 0.67 

56 ±15 0.59 

50.3 ±9.4 0.78 



i.e. in what concerns the behaviour of D(0) in the 
infinite volume the analysis of fitting do not allow 
to conclude anything. In order to distinguish between 
a AjV"^ behaviour or that summarized in equation 
(O simulations with larger volumes are needed and/or 
simulations with larger number of configurations; 
unfortunately, both solutions are beyond our current 
computational power. Of course, a theoretical guidance 
on the volume behaviour of these quantities would be 
also welcome. 



nal analysis of the SU(2) gauge propagator [23|, we start 
assuming that the lattice data has a volume de- 

pendence. This functional form describes quite well the 
lattice D(0) both for the SU(2) and SU(3) propagators. 
For SU(3), the data supports the prediction of i-C. 
£'oo(0) = in minimal Landau gauge, and the Gribov- 
Zwanziger confinement scenario. This is a major differ- 
ence between the SU(2) and SU(3) simulations. What is 
the origin for such difference is an open question which 
deserves further investigations. We would like to recall 
the reader that, in this work, a stability analysis was 
performed and the results look rather robust against a 
change in the volumes used in the fits. Furthermore, the 
sca ling analysis is consistent with the results reported 
in [iSj, where the authors have shown that the infrared 
gluon propagator is compatible with a pure power law 
(p^)^'', with K ~ 0.53. Recall that this previous result 
implies, again, a -Doo(O) = 0. 

Besides the scaling analysis assuming a be- 
haviour, more general ansatze for the dependence with 
the lattice volume were considered. The results of using 
these ansatze are not conclusive but open the possibility 
of having a Doo(O) ^ 0. In this sense, the simulations dis- 
cussed in this article do not provide a definitive answer 
about Doo(O). Clearly, a solution to the puzzle require 
SU(3) simulations with larger volumes and a deeper theo- 
retical understanding of the behaviour of I?(0) and M (0) 
with the lattice volume. 

One issue which was not considered neither here nor in 
the SU(2) simulation [2^, is the effect due to the Gribov 
copies [2^, [13, HH, nil • However, given that the a values 
reported here are several cr's above the critical value of 
0.5 for (Af(0)) and 1 for D{^)lV and (M(0)2), it seems 
unlikely that Gribov copies would change our conclusion. 
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Besides the lattices reported in table U] we have per- 
formed similar analysis on the following two sets (with 
13 = 6.0): i) 16^ x T for T € {16,32,64,128,256}; ii) 
X 32 for L G {16, 20, 24, 28, 32}. The analysis of these 
two sets also gives support for D{0) — 0. Indeed, for the 
first set we have D{0) ~ V'-°°''2'"> (x^/d.o.f. = 0.11) 
and {M(Of) ~ y-i-090(ii) (x^/d.o.f. = 0.27). For set ii) 
the scaling analysis gives D(0) ~ V'"" ''^'''^^' (x^/d.o.f. = 
0.67) and {M{Of) ~ [x^/d.o.f. = 1.68). 



